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Abstract 

The search for an engine cycle with high efficiency, multi-sources of energy and less pollution has led to reconsideration of the Stirling 
cycle. Several engine prototypes were designed but their performances remain relatively weak when compared with other types of 
combustion engines. In order to increase their performances and analyze their operations, a numerical simulation model taking into 
account thermal losses has been developed and used, in this paper, to optimize the engine performance. This model has been tested using 
the experimental data obtained from the General Motor GPU-3 Stirling engine prototype. A good correlation between experimental data 
and model prediction has been found. The model has also been used to investigate the influence of geometrical and physical parameters 
on the Stirling engine performance and to determine the optimal parameters for an acceptable operational gas pressure. 

© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 

The urgent need to preserve fossil fuels and use renew¬ 
able energies has led to the use of Stirling engines, which 
have an excellent theoretical efficiency, equivalent to the 
Carnot one. They can use any source of energy (combus¬ 
tion energy, solar energy...) and they are less polluting 
than the traditional engines. 

Several prototypes were produced, Fig. 1, but their 
outputs remain very weak compared to the excellent 
theoretical yield, [1-5]. In fact, these engines have 
extremely complex phenomena related to the compressible 
fluid mechanics, thermodynamics, and heat transfer. An 
accurate description and understanding of these highly 
non-stationary phenomena is necessary so that different 
engine losses, optimal performance and design parameters 
can be determined. 

Many investigators have studied the effect of some heat 
losses and irreversibilities on the engine performance 
indices. However, they have not calculated the optimal 
performance and design parameters for maximum power 
and efficiency. Popescu et al. [6] show that the most 
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significant reduction in performance is due to the non- 
adiabatic regenerator. Kaushik, Wu and co-workers [7,8] 
have found that heat conductance between the engine and 
the reservoirs, the imperfect regenerator coefficient and the 
rates of the two regenerating processes are the important 
factors affecting the performance of a Stirling engine. 
Kongtragool and Wongwises [9] investigated the effect of 
regenerator effectiveness and dead volume on the engine 
network, heat input and efficiency by using a theoretical 
investigation on the thermodynamic analysis of a Stirling 
engine. Costea et al. [10] studied the effect of irreversibility 
on solar Stirling engine cycle performance; they included 
the effects of incomplete heat regeneration, internal and 
external irreversibility of the cycle as pressure losses due to 
fluid friction internal to the engine and mechanical friction 
between the moving parts. Cun-quan et al. [11] have 
established a dynamic simulation of an one-stage Oxford 
split-Stirling cryocooler. The regenerator inefficiency loss, 
the solid conduction loss, the shuttle loss, the pump loss 
and radiation loss are integrated into the mathematical 
model. The regenerator inefficiency loss and solid conduc¬ 
tion loss are the most important. An acceptable agreement 
between experiment and simulation has been achieved. 
Cinar et al. [12] manufactured a beta-type Stirling engine 
operating at atmospheric pressure. The engine test 
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Nomenclature 

Subscripts 

A 

area, m 2 

c 

compression space 

c? 

specific heat at constant pressure, Jkg -1 K -1 

diss 

dissipation 

c 

^pr 

heat capacity of each cell matrix, WK -1 

d 

expansion space 

£ 

regenerator efficiency 

E 

entered 

M 

mass of working gas in the engine, kg 

f 

cooler 

m 

mass flow rate, kgs -1 

h 

heater 

m 

mass of gas in different component, kg 

irr 

irreversible 

P 

pressure, Pa 

P 

loss 

Q 

heat, J 

Pa 

wall 

Q 

power, W 

r 

regenerator 

R 

gas constant, JkgK -1 

r i 

regenerator cell 1 

T 

temperature, K 

U 

regenerator cell 1 

U 

convection heat transfer coefficient, W m -2 K -1 

S 

outlet 

V 

volume, m 

shtl 

shuttle 

W 

work, J 

T 

total 


indicated that the engine speed, engine torque and power 
output increase proportionally with a rise in the hot source 
temperature. 

Walker [2] mentions other losses but without introducing 
them in the models: the conduction losses in the 
exchangers, the load losses, the shuttle losses and the gas 
spring hysteresis losses. Furthermore, these losses are not 
usually studied in literature because of their complexity. 
Urieli and Berchowitz [13] developed an adiabatic model 
and a quasi-stationary model where they introduced only 
the pressure drops into the exchangers. The results 



obtained by this model are better than those of the other 
models, but remain different from the experimental results. 
Hence, the Stirling engine performance depends on 
geometrical and physical parameters of the engine and on 
the working fluid gas properties such as regenerator 
efficiency and porosity, dead volume, swept volume, 
temperature of sources, pressure drop losses, shuttle losses, 
etc. 

A dynamic model taking into account the different losses 
is developed by the authors and tested using the General 
Motor GPU-3 Stirling engine data, [14]. The results 
obtained proved better than those obtained by other 
models and correlate more closely with experimental data. 
The model is used to determine the losses in different 
engine compartments and to calculate the geometrical and 
physical parameters corresponding to minimal losses [15]. 
An optimization based on this model is presented in this 
article. It will help study the influence of geometrical and 
physical parameters on the prototype performance of a 
Stirling engine and therefore determine their optimal 
values. 

2. Dynamic model including losses 

A second-order adiabatic model has been initially 
developed. The estimated values of the engine parameters 
are obtained and for the sake of validation, the results are 
compared with Berchowitz results under analogous condi¬ 
tions [14]. Afterward, a dynamic model, which takes into 
account the losses in the different engine elements, was 
developed. 

The losses considered in this model are the energy 
dissipation by pressure drops in heat exchangers, energy 
lost due to internal conduction through the exchangers, 
energy lost due to external conduction in the regenerator, 
energy lost due to the shuttle effect in the displacer and 
energy lost due to gas spring hysteresis in the compression 
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and expansion spaces [14]. The mechanical friction between 
the moving parts is not considered. 

The schematic model and the temperature distribution in 
the various engine components are shown in Fig. 2. The 
dynamic model of the developed Stirling engine is based on 
the following assumptions: 


V-h - 


3T r 2 — Tri 
2 


( 5 ) 


According to the flow direction of the fluid, the interface’s 
temperatures are defined as follows [14]: 


TV f = T c if m c _f>0, otherwise TVf — Tf. 


The gas temperature in the different engine elements is 
variable. 

The cooler and the heater walls are maintained 
isothermal at temperatures T pa f and T pah . 

The gas temperature in the different components is 
calculated using the perfect gas law. 

The regenerator is divided into two cells rl and r2, each 
cell has been associated with its respective mixed mean 
gas temperature T r{ and T r2 expressed as follows: 


TV r = V if Wf_ r >0, otherwise Tf_ r = T r _ r . 

7V h = 7V r ifra r _h>0, otherwise TVh = V. 

TV d — T h if Wh-d>0, otherwise TVd = TV 

The regenerator matrix temperatures are therefore 
given by 

depart _ _ 8 (2rl 

d t C pr dF K ’ 


Pi 1 Vrl 

Rm r \ 


( 1 ) 


Pr2 Vr2 

Rm r2 


( 2 ) 


An extrapolated linear curve is drawn through temperature 
values TVi and T t2 , defining the regenerator interface 
temperature TV, TV and T r _ h , as follows [15]: 


T f- r 


3 Vi — T r 2 
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± r— r 


T rl + T r 2 
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(3) 
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Fig. 2. Schematic model of the engine and various temperature 
distributions. 


dT par 2 _ _ 8Q r 2 

d t C pr dt' { ’ 

Taking into account the losses by internal conduction in 

the exchangers: 50 pcdf , 8 2 pc dri> §2 Pc dr 2 > 8 2 pc dh and 
external conduction in the regenerator, the power ex¬ 
changed in the different heat exchangers are given by 

6 Qf = UfApAT^s - Tf) - 5e pcdf , (8) 

S&, = eU rl A parl (Tp,„ - T n ) - 8 ^1, (9) 

5 Q i2 = eUrtA^ATpvi ~ T, 2 ) - , (10) 


8 (2h = pah(Tpah — V) - S() pcdh , (H) 

where s is the regenerator efficiency. 

The heat transfer coefficient of exchanges U { , U r] , U r2 
and U h are available only empirically [13]. The total 
exchanged heat is 

s Q = sg f + 5g rl + 5g r2 + 8g h - 50 pshtl , (12) 

where d0 pshtl is the shuttle loss in the displacer. Consider¬ 
ing the loss due to gas spring hysteresis in the compression 
and expansion space: 8W mc /dt and bJFVd/dt, evaluated in 
[14], the work done in a cycle is 

8W_ dV dV _ VfVe _ 514Vd 
d t c d t d dt d t dt 

The total engine volume is 

V = V c + V f + V rl + F r2 + V h + F d . (14) 

Since there is a variable pressure distribution throughout 
the engine, the compression space pressure P c has been 
arbitrarily chosen as the baseline pressure. At each 
increment of the solution, P c is evaluated from the relevant 
differential equation and the pressure distribution is 
determined with respect to P c . 

The other variables of the dynamic model are given by 
the energy and mass conservation equation applied 
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to a generalized cell: 

s/i i t t* t> dV , ^ d(mT) 

oQ + C p r E m E — CpTsWs = F — + CV——, 


where §(2dissT 8(9 d j ss f T 8(2dissrl T 8(?dissr2 T 8(2dissh’ ^ the 
total heat dissipation generated by pressure drop. 

(15) The mass flow in the different engine components is 
given by the energy conservation equations (17)—(22): 


M — m d + m c + Wf + m r + n%. (16) 

Applying the energy conservation equation to the 
different engine cells and including energy dissipation by 
pressure drop in the exchangers, <5<2 diss , and the other losses 
yields 

— C p r c _fm c s = — + CVU c — ~ 8ITjn-c, 

(17) 


80f — §(?dissf — C p 7>_ r Wfs 


C v Vf d P c 
R d t ’ 

(18) 




+ C p 7>_ r m r , E — C p r r _ r m r is 


Cy V r \ d P c 
R d t ’ 
(19) 


80r2 8^dissr2 T C p T r _ r m r 2 E C p T r _| 1 /71 r 2S 


Cy U r 2 dPc 

_ R dT’ 

( 20 ) 


1 


WcS = 


d '% Kc d ^ + 8 ^" 


i?r c _f \ d t ydtj C P r c _f ’ 


(24) 


1 


Wfs = 


c p Tf_ r 


80f - S<2dissf + c pT c-fhifE: 


Cy Kf dP c 

_ R dT" 


(25) 


_ 1 TxA xb I ^ rp CyV rl dP c 

hlrlS — ^ ( 8(9 r 2 6(2dissrl + C p 2 f- r m r i E 


c p T r _ r 


2? dt 


(26) 


_ 1 u ,S ^ C V F r2 dP c 

h?r2S — t, | 8^/ r 2 8 C^dissr2 “f C p i r—r^r2E 


c p T r _ h 


i? dt 7’ 
(27) 


• _ 1 P/S s/S y rr ■ dm h c y V h dP c \ 

'" hS " c p r h . d ( 5eh 8 2dis S h + c P r r -hnihE dr£ ^ d , )• 

(28) 

with m c s = hif E , mfs = m r | E , m r is = m r 2 E , m r 2 s = «ihE and 
h?hs = w dE . 


8(2h 8(2dissh + CpTj-^m^ CpT^—e/hijS 


CvKhdPe 

i? dt ’ 

( 21 ) 


CpTh-d^d - 8Q ps htl = ( CpPd_ d/+ ~ 8 ITirrd■ 

( 22 ) 

Summing Eqs. (22)-(27), the pressure variation was 
obtained: 


cLPc 

dt 


1 

C v F t 


8 IT' 

7^(80-82 dissT )-C p — 


(23) 


3. Prototype specifications 

The developed model has been tested using data from the 
Stirling engine GPU-3 manufactured by General Motor; 
this engine has a rhombic motion transmission system, Fig. 
1. The geometrical parameters of this engine are given in 
Table 1. The operating conditions are as follows: working 
gas helium at a mean pressure of 4.13 MPa, frequency 
41.72 Hz, hot space temperature T pah = 977 K and cold 
space temperature T paf = 288 K. The measured power 
output was 3958 W, at a thermal efficiency of 35%. The 
independent differential equations, obtained in paragraph 


Table 1 

Geometric parameter values of the GPU-3 Stirling engine 


Parameters 

Values 

Parameters 

Values 

Clearance volumes 


Cooler 


Compression space 

28.68 cm 3 

Tubes number/cylinder 

312 

Expansion space 

30.52 cm 3 

Interns tube 

1.08 mm 

Swept volumes 


Diameter 

46.1 mm 

Compression space 

113.14cm 3 

Length of the tube 

13.8 cm 3 

Expansion space 

120.82 cm 3 

Void volume 


Exchanger piston conductivity 

15 W/mK 

Regenerator 


Exchanger piston stroke 

46 mm 

Diameter 

22.6 mm 



Length 

22.6 mm 

Heater 


Wire diameter 

40 pm 

Tubes number 

40 

Porosity 

0.697 

Tube inside diameter 

3.02nun 

Units numbers/cylinder 

8 

Tube length 

245.3 mm 

Thermal conductivity 

15 W/mK 

Void volume 

70.88 cm 3 

Void volume 

50.55 cm 3 
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2, are solved simultaneously for the variables: P c , m c , T rX , 
W, etc. 

The vector Y denotes the unknown functions. For 
example, Y pc is the system gas pressure in the compression 
space. The initial conditions to be satisfied are noted: 

Y(to) = Y 0 . 

The corresponding set of differential equations is expressed 
as 


d Y 
d t 


F(t, Y). 


The objective is to find the unknown function Y(t) which 
satisfies both the differential equations and the initial 
conditions. The system of equations is solved numerically 
using the classical fourth-order Runge-Kutta method, 
cycle after cycle until periodic conditions are reached. 

To validate the numerical method used in the computa¬ 
tion, the results are compared with those obtained by Urieli 
and Berchowitz [13] under the same conditions (adiabatic 
models) of the GPU-3 engine data. The comparison shows 
a good agreement [14]. 


4. Results of the dynamic model with losses 

It should be noted that all losses have been included in 
the model; the heat flow rate for each component versus 
crank angle is illustrated in Fig. 3. The corresponding 
average power of the engine is equal to 4.27 kW. The 
average heat flow generated by the heater is equal to 
10.8 kW; it yields an engine efficiency of 39.5. The power 
and the efficiency calculated by the model are very close to 
the power and the actual efficiency of the prototype given 
in paragraph 3. 


xIO 5 



The heat flow loss by internal conduction, the energy 
dissipation by pressure drop through the heat exchangers 
and the shuttle heat loss in the displacer are given in Fig. 4. 
The energy lost due to internal conduction is negligible in 
the heater and in the cooler and is about 8.5 kW in the 
regenerator, which represents 35% of the total energy loss. 
This is due to the lengthwise temperature variation, which 
is very significant in the regenerator. The energy lost due to 
dissipation is mainly observed in the regenerator; it reaches 
a maximum of 3.9 kW, with an average of 935 W. In the 
heater and in the cooler it is equal to 26.6 and 123 W, 
respectively. The average heat flow value lost by the shuttle 



Internal conduction loss 

Energy dissipation by pressure drop 

Shuttle heat loss 


Fig. 4. Lost heat flow in the engine. 


xIO 4 



Fig. 5. Lost heat flow by external conduction in the regenerator. 
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effect is about 3.1kW, which represents 13% of the total 
energy loss. 

The energy lost due to external conduction in the 
regenerator is 27 kW, which represents 47% of the total 
losses (Fig. 5), and is very significant and depends mainly 
on the regenerator efficiency. The energy lost due to 
irreversibility in the compression and expansion spaces is 
very low [15]. 

5. Optimization of the Stirling engine performance 

The energy losses occur mainly in the regenerator. They 
are primarily due to the losses by external and internal 


conduction and pressure drop through the heat exchangers. 
The energy lost due to the shuttle effect in the exchanger 
piston is also significant; it is about 13%. The other losses 
are very weak [15]. 

Reduction of these losses can improve the engine 
performance. Such losses depend mainly on the matrix 
conductivity of the regenerator, its porosity, the inlet 
temperature variation, the working gas mass, the regen¬ 
erator volume and the geometrical characteristics of the 
displacer. To investigate the influence of these parameters 
on the prototype performance, we have varied the studied 
parameter in the model each time and have kept the others 
unchanged, equal to the prototype parameters. 



Fig. 6. Effect of regenerator thermal conductivity on performance. 



Fig. 7. Effect of regenerator heat capacity on performance. 
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5.1. Effect of the regenerator matrix conductivity and heat 
capacity 

The performance of the engine depends on the con¬ 
ductivity and heat capacity of the material constituting 
the regenerator matrix. Fig. 6 shows that with an increase 
of the matrix regenerator thermal conductivity leads 
to a reduction of the performance due to an increase 
of internal conduction losses in the regenerator, [14,15]. 
Fig. 7 shows that the engine performance improves 
when the heat capacity of the regenerator matrix 
increases. 

The matrix of the regenerator can be made from several 
materials. The performances of the engine are given 
according to the material type in Table 2. The performance 
of the engine depends on the regenerator matrix material. 
To increase heat exchange of the regenerator and to reduce 
the internal losses by conductivity, a material with high 
heat capacity and low conductivity must be chosen. 
Stainless steel and ordinary steel are the most suitable 
materials to prepare the regenerator matrix. 


5.2. Effect of the regenerator porosity 

The porosity of the regenerator is an important 
parameter for engine performance. It affects the hydraulic 
diameter, dead volume, velocity of the gas, regenerator 
heat transfer surface and regenerator effectiveness; and 
thus affects the losses due to external and internal 
conduction and the dissipation by pressure drop [15]. 

Engine performance decreases when the porosity in¬ 
creases due to an increase in the external conduction losses 
and a reduction of the exchanged energy between the gas 
and the regenerator, Q r (Fig. 8). A low porosity will give a 
better result. 

5.3. Effect of the regenerator temperature gradient 

( Tf- r -T r _ h ) 

Although engine losses increase when the temperature 
gradient of the regenerator increases [15], the performance 
of the engine also increases (Fig. 9). In this case, the 
performance enhancement is due to an increase of the 


Table 2 

Effect of nature of the regenerator metal on the engine performance 


Regenerator matrix 
metal 

Volumetric capacity 
heat (J/m 3 K) 

Conductivity 
(W/rn K) 

Engine power 
(W) 

Engine effectiveness 
(%) 

Exchanged energy in 
the regenerator (J) 

Steel 

3.8465 x 10 6 

46 

4258 

38.84 

441.25 

Stainless steel 

3.545 x 10 6 

15 

4273 

39.29 

448.72 

Copper 

3.3972 x 10 6 

389 

- 

- 

- 

Brass 

3.145 x 10 6 

100 

4080 

34.6 

415.67 

Aluminum 

2.322 x 10 6 

200 

3812 

29.16 

378.03 

Granite 

2.262 x 10 6 

2.5 

4091 

34.51 

430.75 

Glass 

2.1252 x 10 6 

1.2 

4062 

33.85 

427.8 
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exchanged energy between the matrix and the working fluid 
of the regenerator. 

5.4. Effect of the regenerator volume 

To vary the regenerator volume, the diameter is fixed 
and the length is varied or vice versa. When the regenerator 
diameter is fixed at 0.0226 m, the length affects the 
performance. Although the energy exchanged increases, 
the engine power and efficiency reach a maximum. At a 
length equal to 0.01 m, then, the power decreases quickly, 
as shown in Fig. 10. This can be explained by an increase of 
the dead volume. 


When the regenerator length is constant and 
L = 0.022 m, the performance decreases when the regen¬ 
erator diameter increases (Fig. 11). The dead volume and 
the exchanged energy in the regenerator also decrease. 

5.5. Effect of the fluid mass 

An increase of the total mass of gas in the engine leads to 
an increase in the density, mass flow, gas velocity, load and 
function pressure. Therefore, the increase in the total mass 
of gas in the engine leads to more energy loss by pressure 
drop [15]; however, the engine power increases and the 
efficiency reaches a maximum of about 40% when the mass 
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is equal to 0.8 g (Fig. 12). When the mass increases further, 
the decrease of the efficiency is due to an increase of the 
pressure loss and the limitation of the heat exchange 
capacity in the regenerator and the heater. The use of a 
mass of gas equal to 1.5 g in the engine leads to an 
acceptable output and a higher power than in the 
prototype. 

5.6. Effect of expansion volume and exchanger piston 
conductivity 

The expansion volume and the exchanger piston 
conductivity considerably affect the losses due to the 


shuttle effect, which represents 13% of the engine total 
losses. To vary the expansion volume, the stroke constant 
can be maintained and the piston surface can be varied or 
vice versa. When the piston stroke is constant and equal to 
the prototype value of 0.046 m, the effect of the piston 
surface on the performance is given in Fig. 13. When the 
section increases, the engine power increases, but the 
efficiency reaches a maximum. If the exchanger piston area 
is equal to 0.0045 m 2 , a power higher than 5kW and an 
output slightly lower than that of the prototype can be 
reached. When the exchanger piston area is constant and 
equal to the prototype value of 0.0038 m", the effect of the 
stroke variation on the performance is given in Fig. 14. 




220 


130 


140 


100 


60 


20 

-3 


e 


Fig. 12. Effect of fluid mass on performance and engine mean pressure 
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Fig. 13. Effect of exchanger piston area on engine performance. 
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Fig. 14. Effect of exchanger piston stroke on engine performance. 


When the stroke increases, the engine power decreases, but 
the efficiency reaches a maximum. The optimal perfor¬ 
mances are superior to that of the prototype. They are 
obtained when the area and the stroke are, respectively, 
equal to 3.8 10“ m and 0.042m, which correspond to a 
power of 4500 W and an efficiency of 41%. 

The thermal conductivity of the exchanging piston 
considerably affects the engine performances (Fig. 15). A 
weak conductivity reduces the losses due to the shuttle 
effect and consequently increases the engine power and the 
efficiency. 


6. Conclusion 

The Stirling engine prototypes designed have low 
outputs because of the considerable losses in the regen¬ 
erator and the exchanger piston. These losses are primarily 
due to external and internal conduction, pressure drops in 
the regenerator and shuttle effect in the exchanger piston, 
which depend on the geometrical and physical parameters 
of the prototype design. An optimization of these para¬ 
meters has been carried out using the GPU-3 engine data, 
and has led to a reduction of the losses and to a notable 
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Fig. 15. Effect of displacer thermal conductivity on performance. 


improvement in the engine performance. The parameters of 
this prototype were first applied on the developed model; 
the results were very close to the experimental data. Then, 
the influence of each geometrical and physical parameter 
on the engine performance and the exchanged energy of the 
regenerator was studied. 

A reduction of the matrix porosity and conductivity of 
the regenerator increase the performance. A rise of the 
total mass of gas leads to an increase of the engine power 
and function pressure; however, the efficiency reached a 
maximum. When the displacer section increases and the 
piston stroke decreases, the engine power increases, and the 
efficiency attains a maximum. A low conductivity of the 
exchanger piston reduces the losses due to the shuttle effect 
and consequently increases the engine power and efficiency. 

In future, it is hoped to introduce simultaneously in the 
model all the optimal parameters obtained, determine the 
optimal design parameters and consequently the interesting 
performance. 
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